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Evolution of Non-linear Fluctuations in Preheating after Inflation 
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We investigate the evolution of the non-linear long wavelength fluctuations during preheating after inflation. 
By using the separate universe approach, the temporal evolution of the power spectrum of the scalar fields and 
the curvature variable is obtained numerically. We found that the amplitude of the large scale fluctuations is 
suppressed after non-linear evolution during preheating. 
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I. INTRODUCTION 


Preheating after inflation is a crucial stage in the early universe. Fluctuations of matter fields and gravitational fields are 
amplified due to the parametric resonance caused by the coherent oscillation of the inflaton field[l, 2, 3, 4, 5]. One of the 
important question on preheating is how the long wavelength metric fluctuation is amplified by the parametric resonance. The 
linear analysis shows that non-adiabatic modes of long wave fluctuations are amplified and grow during preheating until the non¬ 
linear effect by the second order fluctuations becomes dominant. When the non-linear effect cannot be neglected, the backreaction 
of the fluctuations on the evolution of background quantities becomes dominant and the amplification of long wave modes stops 
[6, 7]. To obtain complete understanding of the evolution of inhomogeneities during preheating, we have to investigate non-linear 
evolution of fluctuations . 

The separate universe approach[8] is an appropriate method to treat the non-linear dynamics of long wavelength fluctuations 
during preheating. This approach neglects fluctuations of which wavelength is smaller than the Hubble horizon scale. For each 
spatial point, the basic equation of the separate universe reduces to that of the Friedmann equation for a homogeneous and isotropic 
flat universe. However, all dynamical variables include non-linear inhomogeneities of which wavelength is larger than the Hubble 
horizon scale. As the separate universe approach includes all order long wavelength non-linear gravitational fluctuations, we can 
apply this method to investigate the backreaction effect on evolution of a Friedmann-Robertson-Walker(FRW) universe [9]. 

Tanaka and Basset[10] investigated preheating using the separate universe approach. They found that initial small fluctuations 
are amplified and evolve to random spatial distribution. In the separate universe, dynamical variables at each spatial point evolve 
independently. Hence, in the early stage of evolution, fluctuations keep their initial spatial distribution and only the amplitude of 
fluctuations grows by the parametric resonance. When the amplitude of the massless held grows to be comparable to that of the 
inflaton field, the non-linear interaction between scalar fields becomes to be dominant and the system enters a chaotic regime. 
Then, the held variables at each spatial point behave as independent random variables. At this stage, the power spectrum of 
fluctuations is same as that of random white noise. 

In this paper, we concentrate on the evolution of the power spectrum of fluctuations during preheating especially on curvature 
fluctuation that gives an impact on formation of large scale structure. By performing the numerical simulation of preheating 
based on the separate universe approach, we aim to understand the feature of evolution of long wavelength non-linear fluctuation 
during preheating. 

The plan of paper is as follows. In Sec. 2, we review the separate universe approach. In Sec. 3, we analytically estimate the 
evolution of the power spectrum of scalar fields and a curvature variable. We present our numerical results in Sec. 4 and Sec. 5 
is devoted to summary and conclusion. We use units in which c = h = SttG = 1 throughout the paper. 

II. SHORT REVIEW OF THE SEPARATE UNIVERSE APPROACH 


The separate universe approach[8] is the lowest order spatial gradient expansion (GE)[11] that is obtained by neglecting all 
terms containing the second order spatial derivative terms in the Einstein equation. All dynamical variables keep the spatial 
dependence that represents large scale inhomogeneities whose wavelength is larger than the Hubble horizon scale. In this 


*Electronic address: nambu@gravity.phys.nagoya-u.ac.jp 
^Electronic address: araki@gravity.phys.nagoya-u.ac.jp 


Typeset by REVT^ 



2 


approach, we assume the following form of the metric : 

ds^ = —N'^dt^ + e^°‘dx^, ( 1 ) 

where iV is a lapse function and a is the logarithm of the scale factor of the universe. For the system with two scalar fields 
t/)i (i = 1, 2), the Einstein equation becomes 
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where Eq. (2) is the Hamiltonian constraint, Eq. (3) is the momentum constraint, Eq. (4) and Eq. (5) are evolution equations for 
the scalar factor and scalar fields. 

To understand how the spatial inhomogeneity is included in dynamical variables, it is convenient to use the Hamilton-Jacobi 
method[ll]. By introducing the Hubble function H = <j> 2 ) that is a function of scalar fields, the above set of equations 

can be written 


= +y(hM. 
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The solution of the Hamilton-Jacobi equation (6) contains two constants of integration di{x),d 2 {x). If these constants do 
not have spatial dependence', the Hubble function H = H{(j)i{x), (j> 2 ix); di, ^ 2 ) satisfies the momentum constraint (7). By 
differentiating H with respect to dj, we have other integration of constants c{x), f{x): 

= p-3c(q;) 

dd,- 
0(12 


( 10 ) 

( 11 ) 


By solving these equations with respect to (f>j, we can obtain the scale factor dependence of the scalar fields: 

</>! =(j)i{a + c{x); f{x), di,d 2 ), 

<p2=(l)2{a + c{x)]f{x),di,d2). ( 12 ) 


The gauge invariant variable that reduces to the spatial derivative of the scalar field on the zero curvature gauge is defined by[12] 


Qi — ^j^adi(y.j 


(13) 


^ The constants dj can have spatial dependence if they satisfy didi + did 2 = 0. However, the inhomogeneity caused by di is related to the decaying 
mode of perturbations and sub-dominant. Hence we ignore this type of inhomogeneity here. 
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and these variables are non-linear generalization of Mukhanov’s gauge invariant variable in the linear perturbation. The scalar 
field equation (5) in the zero curvature slice dia = 0 is 


4>j,aa + ^3 -f — 0, 

and by taking the spatial derivative of this equation, it is possible to show that the variable satisfies the 

following evolution equation; 

Qi,aa + ^3 + Qi,a + ^Qi = 0 ) ( 14 ) 

that has the exactly same form as in the linear perturbation case. However, the coefficients of the equations are spatial dependent 
and this equation describes non-linear evolution of large scale inhomogeneities. The non-linear generalization of the gauge 
invariant variable in the linear perturbation that reduces to the spatial derivative of the three curvature on the comoving slice is 
defined by 


NH. 




TZi = dia - ^T^dtptoi = dia - "''' = dia - ^diH, 
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and this quantity can be written using Qp. 
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By substituting the solution (12), the gauge invariant quantity Q becomes 


Q^ 


( ^l.a (Pl,f \ (d,c{x)\ 
V 02 ,a ^ 2 ,/y \dif\x)) 



and the gauge invariant variable TZi becomes 


TZi = dia — diC + 


Ej 4’j,a4’jJ 
E, </>,-a 


If we take the zero curvature slice dia = 0, 


TZi = 


-diC + 


Ej 4‘j,a4‘j,f 

E, 


(15) 
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(17) 
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From this formula, we can observe that the curvature fluctuation caused by c{x) does not evolve. This mode corresponds to the 
adiabatic mode in linear perturbation. The curvature fluctuation caused by f{x) corresponds to the non-adiabatic mode and the 
evolution of the curvature fluctuation is caused by the non-adiabatic mode of scalar fields fluctuations. 


III. EVOLUTION OF INHOMOGENEITIES IN PREHEATING 


In this section, we consider a specific inflationary model with the inflaton field (j) and the massless scalar field x- The potential 
is assumed to be 


V(^, x) = (20) 

For ^ nipi, the inflaton field oscillate coherently about ^ = 0 and the evolution of the x Held gets the effect of the parametric 
resonance. For g^/X = 2 case, the longest wavelength mode fc = 0 is included in a strong resonance band and the fluctuation 
grows exponentially in time with Floquet index g w 0.238. Thus, the super-horizon mode is amplified by the parametric 
resonance. As we are interested in the non-linear evolution of fluctuations, we consider g^/X = 2 case in which the growth rate 
for k = 0 fluctuation is the largest. 
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A. Analytic approximation 

In preheating model defined by (20), the time averaged equation of state of the inflaton fields is same as that of radiation and 
the evolution of the scale factor is given by 


e 



( 21 ) 


where we choose a = 0 at f = 0. By introducing the conformal time rj = J die “ and the conformal variables <j) = 
the evolution equations of scalar fields become 


^ 0 , 


e°‘(t>,x = e“x. 


( 22 ) 


where ' = d/dg. The curvature variable (19) becomes 


TZi 



(j)'di^ + x'd^x + O 



(23) 


The solution of the separate universe can be written as 

^ = HvJix)), x = xiv,fix)), (24) 


where f{x) is a spatially dependent constant of integration. We obtain perturbative solution by assuming that the amplitude of 
X is small: 


(j) = (j)o + (j)i-{ -, x = -. (25) 

Then the equations for each order of the perturbation become 

+ X^l = 0, (26) 

Xi+g^^lxi=0, ^1 = 0, (27) 

^2 + 3A^o^2 = -g'^^oXi, X2 = 0, (28) 

Xa + 5 ^^oX 3 =-2g^^o02Xi- (29) 

The solution of the background inflaton field <j>Q is given by 

^0 = ccn(A^^^cr7; 1/2) (30) 


where cn is the Jacobi elliptic cosine, c is the initial value of (pQ. We approximate the elliptic function as follows[13]: 

C cos(t), PI w [Fo + Fi cos(2r)], (31) 

Oyr 

r = T = 7.416, Fq = 0.457, Fi = 0.4973. 

The equation of xi becomes the Mathieu equation and xi grows by the parametric resonance: 

X1 = /e''^ cos r, /r « 0.12y (32) 

where / is a spatially dependent constant of integration which specifies the initial amplitude of x- The solution up to the third 
order becomes 


P = 
X = 



/e'^" - C2 



cos(t) + (sin part + higher frequency part), 


c 


fg^^ 





cos(r) + (sin part + higher frequency part), 
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(33) 

(34) 



5 


where ci, C 2 are 0(1) numerical factors. At the time given by 


the amplitude of the higher order perturbation terms grow to be comparable to that of the lower order quantities. At this time, 
fluctuations of x field changes the evolution of (p field through back reaction effect and the exponential growth of fluctuation of 
X-field is shut off. After this time, non-linearities become dominant. 

As the Eq. (22) is a chaotic system, the behavior of the scalar fields becomes chaotic after this time. The time scale (35) 
depends on the initial value f{x) of x-field that is the spatial dependent function. The deviation of the field of neibouring two 
spatial points grows exponentially for g < rj^,. After 77 > ryn., as the system is chaotic, the dependence of the initial condition is 
completely randomized and the spatial distribution of x becomes that of white noise [ 10 ]. 

To evaluate the evolution of the power spectrum of x, we assume the function f{x) that specifies the initial distribution of x 
obeys the random Gaussian statistics: 




ihJL) = fo 


^kik 2 fk = f- 


where kc is a cut off of wave number and fo defines the amplitude of the power at fc = kc- The actual initial power spectrum 
for the x-field expected to be produced during inflation in this model is investigated by Zibin et al.[7] Owing to their result, the 
initial spectrum of x-field for large scales is given by n = —3 and scale invariant. In this paper, as we are interested in feature of 
non-linear evolutions of fluctuations, we does not fix the value of power index of x-field. The power spectrum of x field is given 
by 






where Ci,C 2 are 0(1) numerical factors. At p ~ ry*, the amplitude of the fluctuation with the wave number 

fk^\ f Xc'^\^ ( 2nX^/^cq^\ 

Uj~W 

becomes 0(1) and non-linear. Once the fluctuation with the wave number k^, becomes non-linear, as the system is chaotic, the 
coherence beyond the scale l/A:* is lost and the fluctuation with the wave number fc < A:* is randomized and behaves as white 
noise. Thus, the power spectrum for fc < /c* becomes that of white noise oc k^. The evolution of the shape of the power spectrum 
depends on the value of the initial power index n. 

For n < —3, as the power of smaller k mode is larger, the smaller k mode (large scale) goes to be non-linear first. At 77 ~ rj^:, 
fluctuation with the wave number k < k^, reaches non-linear and the power spectrum becomes ^ k^- At this time, fluctuation 
with the wave number k > k^: stays linear and the spectrum keeps its initial shape ~ Thus, the shape of the power 

spectrum at 77 ~ 77 * is given by 


r^{k) 


{k < Ac*) 
(k > k^f) 


For n > —3, large k modes (small scale) goes to non-linear first. Once the smallest scale of the system becomes non-linear at 
77 ~ 77 *, the coherence of larger scale fluctuation is lost and the spectrum for all wave number changes to be that of the white 


'Px(fc) ~ {k < k^) 


The power spectrum of the curvature variable TZi is given by 


f k\ / k 


r4 


where C 3 , C 4 are 0(1) numerical factors. At 77 ~ 77 *, the power spectrum for n < —3 is 


Vn^ik) 


kf {k < fc*) 

(fc > fc,) 
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For n > —3, the spectrum at 77 ~ p* is 


'PiZi{k) ^ {k <k^f) 

We summarize the result of the analytic estimation of the power spectrum in TABLE I. 


n < —3 


v<v* 

v>v* 

Px 

kr^+3 

{k < k,), (k > K) 

P'R-i 


k^ {k < fe*), fc"+® (fc > fc.) 


n > —3 


V<V* 

v>v* 

Px 

k’^+s 

k^ 

Pr-i 

k«+5 

k^ 


(43) 


TABLE I: Power spectrums of x and TZi for 3-dimensional space. The initial x is scale invariant for n = —3. 

To compare the numerical calculation performed in the following sections, we present the estimation of the power spectrum 
for 1-dimensional case (TABLE II). The initial spectrums are assumed to be 

Vx(xk^+\ Vni0^k^+^. (44) 


n < —1 


77 < 77* 

v>v* 

Px 


k{k< fe.), fc"+i (k > fc*) 

Plti 

k^+3 

k{k< fc.), fc "+3 (fc > fc*) 


n > —1 


7? < 77 * 

77 > 77 * 

Px 


k 

P-Ri 

k«+3 

k 


TABLE II: Power spectrums of x and TZi for 1-dimensional space. The initial x is scale invariant for n = — 1. 


B. Numerical simulation of preheating 


To confirm the analytic estimation of the evolution of the power spectrums, we perform 1-dimensional lattice simulation of 
the separate universe. 

By introducing the momentum variables = (j)/N, = x/N, our basic equations of the separate universe are 


3H^ = ^{P^ + P^) + V{<P,x), 

(45) 

diH = -^{P^di(j) + P^dix), 

(46) 

Gt J-J. 1 y 0 0 \ 

N ^ N ^ 

(47) 

i Pch dV 

(48) 

A - p E - _3i7P _ E 

N~ N~ ^ dx 

(49) 


As the initial condition, we must prepare variables {H, P^, P^, (f),x) that satisfy the Hamiltonian constraint (45) and the momentum 
constraint (46). We adopt the following form of the initial condition : 

H = ho, Px = 0, X = Xo{x), (j) = (j)o, (50) 

where Hq and (j>o are spatially homogeneous constants and Xo(a^) is the spatially dependent function that specifies the initial 
inhomogeneity of the x-field. We prepare xo(a^) as a random Gaussian field: 

Xo{x) = E AkCos{kx + 9k) (51) 

fe >0 


where 9k is a uniform random number in [ 0 , 27r] and Ak is a random number whose probability distribution is given by the 
Rayleigh distribution 


P{Ak) (X Ak exp 


Xo(fc/fcc)” 


(52) 
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where n is the power index of the field. We choose kc = irNg^d/L where L is the size of the calculating region and Ngn^ is 
number of grids. This value is the smallest wave number that corresponds to the size of the system. The value of the momentum 
of the inflaton field is determined by the Hamiltonian constraint (45): 

P^=[6hl-2V{<l>o,Xo{x))f^. (53) 

Once the initial values that satisfy the Hamiltonian constraint and the momentum constraint, these constraints are preserved 
during temporal evolution. 

We then must specify the time slicing condition to evolve the system. The simplest one is the synchronous slicing = 1. In 
this slice, the evolution equations become 

a = H, i7=-i(p| + p 2 )^ (54) 

^ = P^, P^ = -mp^ - ( 55 ) 

X = P^, P^ = -2,HP^-g^\. (56) 

One of the other choice is the comoving slice. In this slice, the Hubble parameter H is treated as a time parameter (t = H) and 
the evolution equations become 


— — = H 
NdH 


NdH 

= p 

NdH 


1 

N ~~ 
1 dP^ 
N dH 

N dH 


= -3HP^ - 
■ = -iHP^ - 


\4>^-g^dx\ 

9^(ty^x- 


(57) 

(58) 

(59) 


The second equation of Eq. (57) determines the lapse function N. In our numerical calculation, we adopt the synchronous slicing 
N = 1 for the sake of the technical simplicity. We also assume that the a is homogeneous at initial time and we have only 
non-adiabatic mode of curvature fluctuation initially. 

We solve Eqs (54)-(56) using the 4th order Runge-Kutta method. 


IV. NUMERICAL RESULTS 


The spatial grid number used is 2^"^ = 16384 and the power spectrums are obtained by taking ensemble average of 100 
different calculations for different initial distribution of x field with the same spectral index. We used the parameter A = 
0.1, g^ = 0.2, do = 0.5, xo = 10■^ ho = 0.03. 


A. flat space case 

We first performed the numerical calculation of the scalar fields system in flat space (22) to investigate the evolution of 
non-linear fluctuation without gravitational effect. We obtained the evolution of the power spectrum of x-held numerically. 

The local evolution of the inflaton field and the massless field is shown in Eig. 1. The spatial derivative of the fields is calculated 
for two different initial values xo = 0-5 x 10“®, 10“^. 
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T] T] 


FIG. 1: The left panel shows the evolution of the scalar fields at a specific spatial point. The right panel shows the evolution of the spatial 
derivative <j)^x and x^x- 


In this calculation, the time scale given by Eq. (35) at which non-linear effect becomes significant, is 


77 * ~ 280. 

For r] < rjt, the amplitude of x field grows exponentially in time by the parametric amplification while the amplitude of (f) field 
stays constant. In accordance with this behavior, the spatial derivative of x and (t> grows exponentially in time. The growth rate 
of (j)^x is greater than that of x^^ and at 77 ~ 77 *, ~ X,x- After 77 ~ 77 *, the system enters non-linear regime and the grow of 

X,^,x,X,x stops. 

We next investigate the evolution of the power spectrum of x for different initial power indices. 

a. n = —3 case Fig. 2 is the evolution of spatial distribution of x- Initial smooth distribution evolves to random one. 
Up to 77 *, the amplitude of x grows keeping its initial distribution. After 77 *, x begins chaotic behavior and the correlation of 
the spatial pattern is lost. Fig. 3 shows the evolution of power spectrum of x. In this case, larger scale has greater power of 
fluctuation. At 77 = 77 * ~ 280, the fluctuation of the largest scale reaches non-linear and the characteristic scale A:* appears in 
the form of the power spectrum. For fc < fc*, the initial power index —3 changes to be 1 that corresponds to white noise. The 
spectrum for k > keeps its initial slope and a bend at k = k^, appears in the power spectrum. As time goes on, the value of 
A:* becomes larger and the spectrum for k < kc becomes that of white noise V (x k when A:* reaches kc- Thus, although the 
large scale fluctuation reaches 0 ( 1 ) and becomes non-linear temporally, its amplitude decreases to be smaller than unity when 
the fluctuation with smaller scale reaches non-linear. 





FIG. 2: Evolution of the spatial distribution of x field for the initial power index n +1 


2. From the left to the right panel, 77 = 0, 280, 310, 400. 
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J_^^^^^_ 

-3 -2 -1 0 

Ink 


FIG. 3: Evolution of the power spectrum of x fisld for the initial power index n + 1 = —2. The spectrums at time 77 = 
O(blue), lOO(cyan), 200(green), 300(yellow), 400(orange), 500(red) are shown. 


b. n = —1 case In this case, the initial power spectrum of the x field is scale invariant. Fig. 4 shows the evolution of spatial 
distribution of x and Fig. 5 shows the evolution of power spectrum of x- Before 77 *, the spectrum evolves keeping its initial power 
index 0 and all scales goes to non-linear simultaneously at p = 77 *. Then the power index changes to be 1 for all k. As the same 
as 77 = —3 case, the power of the large scale fluctuation is suppressed after the amplitude of the fluctuation reaches non-linear. 



FIG. 4: Evolution of the spatial distribution of x field for the initial power index n-F 1 = 0. From the left to the right panel, rj — 0, 260, 300, 400. 



-3 -2 -1 0 

Infc 

FIG. 5: Evolution of the power spectrum of x fisld for the initial power index n -\- 1 = 0 at 77 

O(blue), lOO(cyan), 200(green), 300(yellow), 400(orange), 500(red). 
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c. n = 0 case In this case, the initial spectrum is same as that of white noise and the small scale reaches non-linear hrst. 
Fig. 6 shows the evolution of spatial distribution of x Fig. 7 shows the evolution of power spectrum of x- As the initial 
spectrum is same as that of white noise, the power index does not change during whole evolution. 



X{x) 



0 50 100 150 200 250 


xi^) 



0 50 100 150 200 250 


FIG. 6: Evolution of the spatial distribution of x field for the initial power index n-F 1 = 1. From the left to the right panel, r/ — 0, 100, 200, 300. 



Ink 

FIG. 7: The power spectrum of x fisld for the initial power index n -\- 1 = 1 at j; 

O(blue), lOO(cyan), 200(green), 300(yellow), 400(orange), 500(red). 


d. n = 1 case In this case, the initial power index is larger than that of the white noise. Fig. 8 shows the evolution of spatial 
distribution of x Fig. 9 shows the evolution of power spectrum of %. After the small scale reaches non-linear, the initial 
power index 2 reduced to be 1. 



FIG. 8: Evolution of the spatial distribution of x field for the initial power index n-F 1 = 2. From the left to the right panel, r/ — 0, 260, 300, 400. 
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Ink 

FIG. 9: Evolution of the power spectrum of x fi^ld for the initial power index n + 1 = 2 at 

O(blue), lOO(cyan), 200(green), 300(yellow), 400(orange), 500(red). 


B. Cosmological case 

We next solve Eqs (54), (58), (56), and investigate the evolution of the long wavelength metric fluctuations. The behavior of 
variables at a specific spatial point is shown in Fig. 10 and Fig. 11. 




FIG. 10: In the left panel, evolution of scalar fields at a specific spatial point is shown. In the right panel, evolution of the spatial derivative of 
the scalar fields is shown. 
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FIG. 11: Evolution of the curvature variable TZi at a specific spatial point. 


The characteristic cosmic time when the system becomes non-liner is given by 




1 


2 


1460. 


Until this time, x,x grows exponentially in conformal time and also grows exponentially in accord with the growth of x,x 
with larger growth rate than that of x,x- We must take care that the amplitude of the spatial derivative of the scalar fields does not 
reach unity at f ^ f*. This is due to the suppression of the cosmic expansion but the non-linear regime begins at f*. At f f*, 
X,x <t>,x and the back reaction of x on the evolution of (p begins to be significant and the growth of the fluctuation caused by 
the parametric resonance stops. The evolution of metric variables H, a, TZi is same as the evolution of x,x and they grow until 
t when the back reaction cannot be neglected. 

Evolution of the power spectrums of x and TZi depends on the initial power index of the x field. The evolution of the power 
spectrum of x is qualitatively same as the flat space case. For n = — 3 (Fig. 12), until f ~ f*, the spectrum of TZi keeps its initial 
flat shape At f ^ f*, the large scale fluctuation of x field goes to be non-linear and the bend of the spectrum of TZi at k = fc* 
appears. The spectrum of TZi becomes k for k < k^,. The amplitude of the curvature variable TZi stays less than unity during 
whole evolution. 




In/c \iik 


FIG. 12: Evolution of the power spectrums for the initial power index n-\-l = —2 case. The left panel is 'P^ and the right panel is . Time 
is f = O(blue), 400(cyan), 800(green), 1200(yellow), 1600(orange), 2000(red). 


For n > —1 (Fig. 13, Fig. 14, Fig. 15), until t ^ f*, the spectrum of TZi keeps its initial shape At f ~ f*, the small scale 
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fluctuation of x field goes to be non-linear and after this time, the spectrum of TZi becomes k for all wave number. The amplitude 
of the curvature variable TZi stays less than unity during whole evolution. 




Ink Ink 

FIG. 13: Evolution of the power spectrums for the initial power index n -F 1 = 0 case. The left panel is 'P^ and the right panel is Time 
is f = O(blue), 400(cyan), 800(green), 1200(yellow), 1600(orange), 2000(red). 




Ink In/c 

FIG. 14: Evolution of the power spectrums for the initial power index n -F 1 = 1 case. The left panel is Vx and the right panel is 'Pui ■ Time 
is f = O(blue), 400(cyan), 800(green), 1200(yellow), 1600(orange), 2000(red). 







14 



FIG. 15: Evolution of the power spectrums for the initial power index n + 1 — 2 case. The left panel is and the right panel is "P-r.. Time 
is f = O(blue), 400(cyan), 800(green), 1200(yellow), 1600(orange), 2000(red). 



V. SUMMARY AND DISCUSSION 


In this paper, we investigated the evolution of long wavelength fluctuations during preheating after inflation. By using the 
separate universe approach, we obtained the evolution of the power spectrum of long wavelength fluctuations numerically. During 
the linear stage of the evolution, the fluctuation of the massless field grows exponentially in time by the parametric amplification, 
but the power spectrum keeps its initial shape. When the fluctuation of the massless field becomes non-linear, the amplification of 
fluctuations stops by the back reaction effect of massless field on the inflaton field. After this time, the system enters the chaotic 
non-linear stage and the shape of the spectrum changes. For large scale mode, the amplitude of the fluctuation is suppressed 
because the statistical nature of the large scale mode becomes that of white noise after the non-linear evolution. The evolution of 
the curvature variable that has an important role in the cosmological scenario is qualitatively same as the evolution of the scalar 
field fluctuations. Hence we do not expect the significant effect of the parametric amplification during preheating on the evolution 
of the large scale fluctuations of metric variables. On the other hand, for small scale fluctuations, the power grows for larger k 
and this leads to the possibility of the formation of the non-linear structures such as primordial black holes[14]. 
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